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Abstract 



At sufficiently high baryon density, a quark matter is expected to become a 
color superconductor because of the pairing forces mediated by gluons. The 
theoretical aspect of this novel phase of the strong interaction is reviewed with 
the emphasis on the perturbation theory of QCD at high chemical potential 
and low temperature. The derivation of the scaling formula of the supercon- 
ducting transition temperature at weak coupling is explained in detail. The 
Ginzburg-Landau theory of the color superconductivity is also discussed. 
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FIG. 1. The phase diagram of QCD. 



I. INTRODUCTION 



The properties of quantum chromodynamics, QCD, at nonzero temperature and nonzero 
baryon density have become a active area of research involving both high energy physicists 
and nuclear physicists. The asymptotic freedom of QCD is expected to release its fun- 
damental degrees of freedom, quarks and gluons, from the color confinement under these 
extreme conditions. New new phases of the strong interaction, the quark-gluon plasma 
(QGP) at high temperature and the color superconductivity (CSC) at high baryon density 
[1] [2] [3] [4] [5] [6] will emerge. The speculated phase diagram of QCD is shown in Fig.l. 
There are two phases of the color superconductivity, 2SC and CFL. The former corresponds 
to Cooper pairing among the two flavors of quarks, u and d, while s quark is too massive to 
participate. The latter occurs at higher chemical potential and involves s in Cooper pairing 
with the ground state color-flavor locked [6]. 

Experimentally, the QGP phase may be reached by colliding two heavy nuclei at suffi- 
ciently high energy. The recent discovery of the jet quenching effect in gold-gold collisions 
versus deuteron-gold collisions at RHIC (Relativistic Heavy Ion Collider) suggests that such 
a novel phase might have been produced post gold-gold collisions [7]. The color super- 
conductivity phase, though not likely to be created at RHIC, is speculated to exist inside 
a compact stellar object, say a neutron star. Its observational aspects, however, is still 
primitive at this time. 

Considerable theoretical endeavors have been expended in past years to explore different 
regions of the QCD phase diagram. Along the temperature axis or slightly off it with a 
small chemical potential, accumulated results of lattice simulation of QCD have reached a 
consensus that a deconfinement phase transition, accompanied by the restoration of chiral 
symmetry, occurs at T ~ 150MeV. But present methodology of lattice simulation become 
inadequate when the chemical potential becomes large because of the fermion sign problem. 
On the other hand, the pairing instability that underlying the color superconducting phase 
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transition extends to an arbitrarily high chemical potential and perturbative techniques are 
available when the chemical potential becomes much higher than Aqcd- First principle 
calculations have been made along this line and the scaling formula relating the transition 
temperature or the energy gap to the chemical potential and the running coupling constant of 
QCD has been derived systematically [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20]. 
Even though the condition for the perturbative treatment may not be implemented in nature, 
these result serves a rigorous proof as to the existence of the color superconducting phase 
of QCD at sufficiently high baryon density. The region of moderate baryon density, which 
represents the most interesting case from the observational prospect, is unfortunately too 
difficult to warrant a first principle investigation with present analytical and numerical 
techniques. One is compelled to resort to various effective actions a la Nambu-Jona-Lasinio, 
regarding the parameters involved phenomenological. 

In this lecture, I will focus on the color superconducting phase of QCD, especially its 
perturbative region. Several salient features of the color superconductivity will be summa- 
rized in the next section. A perturbation theory of CSC at ultra high chemical potential, 
along with the derivation of the scaling formula of the transition temperature, will be re- 
viewed in sections III and IV. A nontrivial application of the perturbation theory to the 
crystalline color superconductivity will be presented in section V. Section VI is devoted to 
the Ginzburg-Landau theory of the color superconductivity, which is less correlated to other 
sections. Some outlooks of the field together with comments on higher order corrections 
and non-standard pairing possibilities are presented in the final section. As only selected 
topics are discussed in this lecture, the interested readers are recommended to consult more 
comprehensive review articles in the literature [21] for the materials not covered here. 



II. THE COLOR SUPERCONDUCTIVITY OF QCD. 

A. The pairing forces in QCD. 

The thermodynamics of a quark matter at nonzero temperature T and nonzero baryon 
density is described by QCD Lagraingian with a Euclidean time < r < j3, f3 = (/c^T) -1 
and a nonzero chemical potential /x 

1 d 

C = — trF* F l — ■?/>( — igA^)il) + fiip^'ip — ipmip + renormalization counter terms, 

2 uXfj, 

(2.1) 

where = A l ^T l , F„ v = §^ - §^ - ig[A^, A u ] with T l the SU(N C ) generator in its 
fundamental representation, and ip is a Dirac spinor with both color and flavor indices. The 
mass matrix m of quarks is diagonal with respect flavor indices. The chiral limit m — > is 
assumed throughout the lecture except in the sections V and VII when exotic pairing states 
are considered. The renormalized coupling constant appropriate for low temperature and 
high density is given by the running coupling constant evaluated at the chemical potential, 
i.e. 
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FIG. 2. The one-gluon exchange vertex. 
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where N f < ^N c and /i >> Aq C d- Although N c = 3 for QCD, most of the formulations 
developed are valid for arbitrary N c . We shall not fix the value of N c unless it is necessary, 
say in the context of the pairing symmetry in the superconducting phase. 

Perturbatively, the diquark interaction is dominated by the process of one- gluon ex- 
change, as is shown in Fig. 2. The amplitude is simply that of the one-photon exchange 
in QED multiplying the group theoretic factor T t 1 T { 2 , which can be decomposed into a 
color anti-symmetric channel (anti-triplet for SU(3)) and a color symmetric one (sextet for 
SU(3)) [3], i.e. 



Nc + l 



N c -1 



(6 



Ci Ci rc„C 2 



+ 5 



c, c 2 rc„ci 



(2.3) 



As both the electric and magnetic parts of the one-photon exchange are repulsive between 
two electrons flying in opposite directions, the one-gluon exchange interaction between two 
quarks flying in opposite directions is attractive in the color anti-symmetric channel because 
of the negative sign of the group theoretic factor (first term on r. h. s. of (2.3)). 

The formula (2.3) can be generalized to the two quarks in arbitrary irreducible represen- 
tations Ri and R 2 of the gauge group and we have 



(T l R2 Y' 2C2 = E G(R\R 1 ,R 2 )V C '^ C1 > C2 (R\R 1 ,R 2 ) 

-RC-Ri®R2 



(2.4) 



with the group theoretic factor 



G(R\Ri, R 2 ) — -jipR — C Rl — Cr 2 ), 



(2.5) 



where Cr is the second Casmir of the representation R and V(R\Ri ® -R2) projects out the 
irreducible representation R out of R± <g> R 2 and satisfies 



(R\Ri, R 2 ) 



V 



C-, ,c,;ci,C2 



(R\Ri, R 2 ) 



for R' = R 
for R' ^ R 



(2.6) 
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The channel with negative G(R\R±, R 2 ) is attractive. The derivation of (2.4) is shown in 
Appendix A. 

Another pairing force, which is expected to dominate at moderate baryon density, stems 
from the instanton-induced interaction among light quarks discovered by 't Hooft [22] . For 
two quark flavors, it is of NJL type and reads [5] 

CcS - = G ( ~ 8N*(N e - 1) [(^T l A ^r 2 T l A C$) + ^Ct 2 T a1 ^t 2 T a15 C^} (2.7) 
+ 16N *( Ne + 1 ^Cr 2 T l s a^)(i}r 2 T l s a^C^, 

where C is the charge conjugation matrix, T l A ^ is proportional to SU(N C ) generator in 
color antisymmetric ( symmetric ) channel and the second Pauli matrix r 2 acts on isospin 
indices of ip. The coupling strength G > can be expressed explicitly in terms of the 
running coupling constant, the chemical potential and the critical chemical potential for 
chiral symmetry restoration [23]. Like one-gluon exchange, the pairing channel is also color 
antisymmetric. 



B. The transition temperature and the energy gap at weak coupling. 

Although the running coupling constant becomes weak for fj, » Aq CD , the infrared 
divergence, because of the long range propagation of magnetic gluons introduces complica- 
tions to the perturbative expansion. Consequently, the scaling formula for the energy gap 
and the critical temperature differs remarkably from that of a BCS superconductor, which 
is driven by a short range pairing force via one-phonon-exchange. Through the efforts of 
several groups, the first-principle formula for the pairing temperature has been derived. It 
reads 

k B T*& = cc'c"c.A 5 e-f[l + O(glng)], (2.8) 
where the coefficient of the non-BCS exponent 



6A" C P . . 

K =fefr (2 - 9) 



was first obtained in [8], the pre-exponential factor 



c = 1024V2n 3 N J (2.10) 



was found in [9] and [11], and the factor 

d = 2e 7 (2.11) 
with 7 the Euler constant was found in [11] and [16]. The factor 

c" = exp[-^(7r 2 + 4)(iV c -l)] (2.12) 
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comes from the non Fermi liquid behavior of the quark self energy and was determined in 
[14] [19] (The existence of this correction was suggested in [8]). The angular momentum 
dependent factor cj was derived in [16] and is given by 

cj = e~ 6sj (2.13) 
for equal helicity pairing with J > and 

cj = e~ 6s 'j (2.14) 
for cross helicity pairing with J > 1 where 

[0 forJ = , n . 

aj = t£Ui forJ>0 ( 2 - 15 ) 

and 

The superconducting transition temperature T c is the highest pairing temperature that 
is consistent with the pairing symmetry. We have T c = T^ ^ for pairing among different 



flavors and T c = T^ ir 1 - ) for pairing within the same flavor [24] [25] 

The gap energy below T c depends on the Matsubara frequency and reads 



at T = and \u\ > A . Its maximum magnitude A for J = is related to T c via 

A J 7re -7 for two flavor pairing 
jfc B T c \ 2~3 7re~ 7 for color- flavor locking, N c = Nf = 3. 



(2.18) 



C. The pairing symmetry. 

The gap parameter in the super phase, depending on the relative momentum, helicity, 
and color-flavor indices of pairing quarks, possesses rich structures under a spatial rotation 
or a transformation of the internal symmetry group 

SU(3) C x SU(3) R x SU(3) L x U(1) B (2.19) 

for three colors and three flavors in the chiral limit. To determine the energetically most 
favored pairing symmetry, we consider first the representation contents of color SU (3) and 
flavor SU(3), consistent with the signature under a simultaneous interchange of the color- 
flavor indices of the pairing quarks. 

According to the angular momentum dependence of the pairing instability shown in the 
previous subsection, the most favorite diquark condensate should carry zero angular momen- 
tum and pair quarks of equal helicity and acquires a negative sign upon interchanging the 
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momenta and spins of the two quarks. Since the condensate should be antisymmetric with 
respect to interchanging all quantum numbers, the sign under a simultaneous interchange 
of their color-flavor indices ought to be positive, i.e. 

< m% >=< nn > ( 2 - 2 °) 

If we define a SU(3) rotation among flavor indices to be the conjugate of that among color 
indices, the representation contents of the diquark condensate are 

(3 C , 3/) 8(60,6/), (2.21) 

which can be decomposed further into a set of irreducible representations under a simulta- 
neous SU(3) rotation of both color and flavor indices. It was found in Refs. [6] [27] that the 
condensate which minimize the free energy implements the unit representations of this diag- 
onal SU(3) rotation. This pairing symmetry is referred to as the color-flavor locking (CFL). 
The diquark condensate of CFL pairing takes the form 

< r/M >= Mmi - mi) + umi + mi)- ^ 

The two invariants 4>a{S) correspond to the two unit representation, one from each term of 
(2.21). The nonlinear coupling between these two unit representations makes (ps 7^ even 
though the pairing force is within the color antisymmetric channel. Both the gap equation 
[6] [28] and the Ginzburg-Landau theory [29] show that the magnitude of (ps is suppressed 
relative to that of <pA by one order of the expansion parameter involved and we shall come 
back to this point in section VI. 

The CFL condensate breaks the symmetry (2.19) of QCD to 

SU(3) C+R+L (3) x Z 2 (2.23) 

with SU{3) c+r+l standing for a simultaneous SU(3) rotation of color indices, right hand 
flavor indices and left hand flavor indices. Like the standard model of electroweak interaction, 
it is possible to factor out a unbroken U(l) group out of the electromagnetic U(l), which is 
a subgroup of the flavor SU(3), and one of U(l) subgroup of the color 577(3) [6] [30]. The 
gauge potential corresponding to this unbroken U(l) plays to role of the electromagnetic 
field in the super phase. This, together with the energy gap, makes the physics of CFL 
condensate similar to that of an insulator regarding the new electromagnetism [31]. 

Since both right hand flavor indices and left hand ones are locked to the color indices. A 
spontaneous chiral symmetry breaking is thereby induced in the CFL condensate with the 
Goldstone bosons represented by the meson octet. The mass spectrum of the mesons in the 
presence of bare quark masses was computed in [32] [33] and they found that the order of 
the meson mass hierarchy are reversed from that in vacuum with rj' the lightest member of 
the family. 

III. THE PERTURBATION THEORY OF ONE GLUON EXCHANGE. 
A. The perturbation theory of the transition temperature 

The perturbative expansion of the pairing temperature refers to the expansion of In -j—^ — 

™B pair 

according to ascending powers of the running coupling constant g, i.e. 



7 



(p V) (p' V') 

"XX" 

(-p-v) (-p'-v) 

C S C.' SL 

2 2 11 

FIG. 3. Proper vertex function, Y C }' C J' C1 ' C2 (u' ,p'\v,p). 

\n—^ = - + \ + \'g + 0(g 2 ) (3.1) 

where the coefficients A, A', may carry powers of In g. The O(^) term represents the 
leading order, the 0(1) term represents the sub- leading order and the 0(g) term represents 
the sub-sub-leading order and so on. It follows from (2.9) that k = /Jgfy and 

A = 5 In g - ln(2048v / 2A r 7 l ) - 7 + — + - - - In Cj (3.2) 

3 16 

The derivation of (2.9) and (3.2) will be the main subject of the rest of this section and that 
of the coming section. 

To generate the perturbation series (3.1) systematically, we start with the proper vertex 
function, T, corresponding to the scattering of two quarks with zero total energy and mo- 
mentum. The Matsubara energies of incoming and outgoing quarks are denoted by ±u and 
±i/, respectively. Similarly, ±p and ±p' label the incoming and outgoing momenta. Each 
of the superscripts q, i = 1,2, denote the color associated with each leg and the subscripts 
s, which label the states above or below the Dirac sea, are either + or — . Primed vari- 
ables are outgoing and unprimed incoming. This labeling scheme is illustrated in Fig. 10. 
Without special declaration, we shall reserve the Greek letters v = 2irkBT(n + 1/2) for the 
discrete Matsubara energy of a fermion and uj = 2nkBTn for that of a boson. The proper 
four-fermion vertex function satisfies a Schwinger-Dyson equation which, with all indices 
suppressed, may be written as, 

T(u',p'\u,p) = f(u',p'\u,p) + W [ ^- 3 K(v',py',q)T(v",q\v,p), (3.3) 

where T represents the two particle irreducible part (2PI) of a vertex diagram. The kernel 
has the explicit form, 

K^z^y^p) = ^:^^:^(^^p>,p)5 <S1 (^)5 4s2 (-p - v \ (3.4) 

where S s i s (p, v) denotes the full quark propagator with momentum p and Matsubara energy 
v. In order to facilitate the partial wave analysis we found it convenient to associate the 
Dirac spinors u(p) and v (p), which satisfy the Dirac equations (74^ — 27 • p)u(p) = and 
(74P — ?7 • p)v(p) = 0, to the quark-gluon vertex instead of to the quark propagator. Thus 
the vertices written in (3.3) are of the form, 
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r\ + ... 

FIG. 4. The Schwinger-Dyson equation. T is represented by single hashed vertices and T is 
represented by double hashed vertices. The full quark propagator is represented by a solid line and 
the bare quark propagator by a thin line. Gluon propagators, including hard thermal loops, are 
represented by curly lines. The first two orders in the expansion of T and the full quark propagator 
are given. 

r^^sa = U 1 {s' l ,p')U 5 {s' 2 , -p')T lS ,a(sUa(si,p)Uf3(s2, -p), (3.5) 

with the vertex function T^ap given by conventional Feynman rules and U(+,p) = u{p) 
or U(—,p) = v(—p), respectively. However, since the quarks are massless, 75W = — u and 
7 5 t> = — v, to simplify notation we may identify U(s,p) = u(sp). The diagram of the Dyson- 
Schwinger equation, the first few diagrams of the Y^ap and S s ' jS (v,p) are displayed in Fig. 
4 

The next step is to decompose T into irreducible representations of SU(N C ) by either 
symmetrization or antisymmetrization among the initial and final color indices. We have, 

+ \ (6< c ^ - 5 c '^5^)rl s ,, su jv f y, \u,p), 

For the pairing instability, we need only to focus on the attractive antisymmetric chan- 

vi ( 

s[,s' 2 ,si,S2 y 



nel. Proceeding with the partial wave analysis, we expand Tf, s , s (i/,p'\u, p) in terms of 



Legendre polynomials ( for equal helicity pairing), 

r£, S2 , sl , s >',P>,p) = £(2J + l)%t S2lS1)S2 (^,P>,p)^(cos^. (3.7) 

j 

with 9 the angle between p and p . Using a similar expression for Tf, s ( y h / ',p'\h',p), we derive 
from (3.3) the Dyson-Schwinger equation satisfied by 7^ a ^. si S2 (z/',p'|z/,p): 

74, S2 , SllS2 (^,P>,P) = ll^s^W ,P>,P) ( 3 - 8 ) 

[ /-00 

+ 13 ^^ 1 , S2 ,4^ S2 '( ^/ ^p> , ^9)74^4^ S1 , S2 (^^gk,p), 

where the kernel Kl , c has the form, 
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K^^VJfrp) = Wl g 2 S s , {si (p, v)S s ,, S2 (p, -v). (3.9) 

The Dyson-Schwinger equation (3.8) is of the Fredholm type, the stability of whose solution 
is governed by the Fredholm determinant [35] 



(3 ^ Jo 



OG 



V = ^{\-K) = \-- I dpK slS2 , slS2 (v,p\is,p) (3.10) 



1 _ poo poo 

+ w dp l dp ' 

v,v ';si,S2;s 1 ,s 2 



(v,p\v,p) K Sl (v,p\i/,pf) 

K s' v s' 2 ; Sl ,s 2 ty, P'W, P) Ks'^s'^ ty, P'W, P') 



In terms of the eigenvalues, E n , of the kernel, defined by the solutions of, 

Ef sltS2 (v,p)=k B T / ^K^ M (uM^M^ 2 W^ (3-11) 

v,s\,s 2 

we have 

V = H(l-E n ). (3.12) 

n 

Since V appears in the denominator of the solution to the Fredholm equation, its roots 
represent the formation of a Cooper pair of quarks at the transition point to the super 
phase. At weak coupling and sufficiently high temperature, all eigenvalues are small (~ g 2 ). 
As the temperature is lowered some of the eigenvalues may cross one because of the Fermi 
sea and the attractive interaction. The pairing temperature at an angular momentum J, 
T^l, is the highest temperature when the corresponding highest eigenvalue crosses one, i.e. 

E mL \ T=T (J) = 1- (3-13) 

pair 

The superconducting transition temperature is the highest one among the pairing temper- 
atures of different angular momentum, i.e. 

T c = max(7£5>,VJ) (3.14) 

For equal helicity pairing, we have T c = T^°\ 

The Fredholm determinant corresponds to the sum of closed diagrams and is formally 
gauge invariant. The scale parameter Aqcd of (2.2), however is gauge dependent. Such a 
dependence should be eliminated order by order upon matching the infrared region and the 
ultraviolet region of the diagrams. Since the uncertainty of A represents an order 0(g 2 ) 
correction to the running coupling constant (2.2), it contributes to the sub-sub-leading term 
of (3.1) and is beyond the scope of the present lecture. 



B. Hard Dense Loop resummation of the gluon propagator. 

A Fermi sea of quarks is highly polarizable in an external gauge field and it in turn 
screens the strength of the external field. This is how the infrared divergence of QCD is 
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regulated at high baryon density and the effect is included in the hard dense loop (HDL) 
resummed gluon propagator, which takes the form V® b u (k,u) = 5 ab V^ u (k,u) with 



X#(k,u;) 



-id 



k 2 + to 2 + a M (k, u) 



fa 



kikj \ 



(3.15) 



45 



ab 



k 2 + a E (k,u) 



(3.16) 



and X>4j (k, uj) = in Coulomb gauge, where icu is the Matsubara energy of the gluon. The 
ultraviolet divergence has been renormalized. For k << /i, uj « //, the magnetic self-energy 
is given by 



a M (k,oj) = m 2 D f M (j) 



and the electric one by 



with 



a E (k,uj) = m 2 D f E (j) 



f M (z) = \ f(l + z 2 ) tan" 1 - - z 
z L z 



(3.17) 



(3.18) 



(3.19) 



= 1 - ^tan^ 1 - < 1 

z 



and 



2 Nfff 



71* JO 



dqq- 



1 



2, ,2 



•pip-n) _|_ x 



2tt 2 



(3.20) 



(3.21) 



the Debye mass. As we shall see in the next subsection, the pairing region of energy and 
momentum, u « k saturates the upper bound of /m(z) and f E (z), and we have effectively 



7T 



a M (k,uj) ~ ~j m D~r~ 



k 



and 



(3.22) 



(3.23) 



The lack of screening in the static limit, uj = gives rise to the forward singularity discussed 
below. 

The HDL resummed gluon propagator in the covariant gauge takes the form 

<*) = 775 ■ "'„, t P^ - ^ ^ 77 ,A " (3.24) 



K 2 + a M (k,uj) K 2 [A; 2 + (T S (A;,^)] "" 

L = 5 

is the gauge parameter 



where = ^ - Uj, = Pjj = PL = 0, - . „„ 



x 2 



(if 2 ) 2 

p£, K 2 = A; 2 + cu 2 and £ 
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C. Forward singularity and the leading order contribution. 



For the one gluon exchange of Fig. 2, incorporating hard dense loops in the gluon 
propagator, we find the s-wave amplitudes of the scattering between two quarks of positive 
energy states 

rt=l + (v',p'\v,p) = ~j(l + j;) /^(cosfl) (3.25) 



x 



3 _ cos £_(l + cos £)jrf! 1 + cos^ 



[u - u'f + \p - pT + m%f M (^) \p- pf + mlf E (^ 



In the static limit, v = u', and p' = p, the first term of the integrand, the color magnetic in- 
teraction becomes singular in the forward direction, and the integral diverges logarithmically 



and the main contribution to the integral comes from the region where 



in the limit — > 0. If this is the dominant pairing force, we expect that v ~ v' ~ ^ c 



\V ~P\ 3 ~ ^™ 2 dW - v \ ( 3 - 26 ) 

The phase space of pairing is restricted to the vicinity of the Fermi surface with \p' — p\ ~ 
k B T c . Anticipating ~ e~s, we have 

Q„ g -%(™±f (3.27) 



and 

(p' -p) 2 



{p 1 — p} 2 

The integration (3.25) can be approximated by 



« 1. (3.28) 



7S° ++ (^v.p) - -4^ (' + jd log ^r (3 ' 29) 

for v' ~ v ~ /c_bT c , where 

V 2 / 256ttV' 

In terms of the bare quark propagator 



Ss>,s(P^)= ■ i6s ''\ , (3-30) 
iv — sp + fi 



the kernel of the Dyson-Schwinger equation reads 
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The contribution from the states below Dirac sea, being far below the Fermi level, will be 
suppressed by two powers of and will not be considered in this lecture. Accordingly, 
the s-subscripts which label the states above or below the Dirac sea will be suppressed in 
the subsequent discussions. 

The quark propagators in the kernel K J=0 makes the integration over p and the sum- 
mation over v logarithmically enhanced for ksT << p. If the pairing force were free from 
the forward singularity, this would be the only source of the logarithm. The powers of the 
logarithm would not grow with the ascending powers of g in the expansion (3.10). The 
scaling behavior ksT c = cexp ( _ p-) is expected. The exponent (leading term) and the pre- 
exponential factor (the sub-leading term) can be fixed by carrying out the expansion (3.10) 
to the second power of the kernel. The forward singularity, however, introduces another log- 
arithm pertaining to the pairing force and its power will grow with that of g. Consequently, 
the non BCS scaling (2.8) is expected. All terms of the expansion (3.10) have to be retained 
and the determination of the exponent and the prefactor of (2.8) becomes non-trivial. 
To make the integral equation (3.8) tractable, the following approximations are made: 
1). A cutoff vq with k B T c « u « p is introduced for the Matsubara energy and the 
vertex function (3.25) is replaced by 

7 J= V,P>,p) * fl + 1) log-A^o - W\)0("o - M), (3.32) 



12pp' V NJ \v - z> 

2) . The logarithm in (3.32) is approximated by 

In — 1 - ~ ln-rj-r (3.33) 
\v ' — v\ |z/>| 

with |z>>| = max(|z>'|, \v\) [8]. 

3) . The summation over discrete Matsubara energies is replaced by an integral over 
continuous Euclidean energy with an infrared cutoff, i.e. 

^(■■0*/, >T^(-) (3 ' 34) 

The integral equation (3.8) under these approximation becomes separable with respect 
to the variable p and can be reduced to that with one variable, u, in terms of the ansatz 

f(u,p) = -(j>(u)e(u -\u\) (3.35) 



p 



and we have 



rb i i 

4>{x) = k 2 / dx'mm(x,x')(f)(x') (3.36) 

J a 



where 



k = 



9' 



(1 + ^), (3.37) 



2An 2 E V N, 
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a — In I, b — In -, and x = In I with 



= (?) 



N f ^ 9 5 k B T 
2 / 256ttV 



and 



2 / 256vrV 



(3.38) 



(3.39) 



Notice that the kernel is symmetric with respect to v and — v and only even functions 
4>{—v) = 4>{v) are considered. The odd functions which vanishes at v — will not contribute 
to pairing since 0(0) = 0. 

Differentiating both sides of (3.36) twice with respect to x, it is reduced to a differential 
equation with trigonometric solutions, i.e. 



d 2 (j) 
dx 2 



+ k 2 



0. 



(3.40) 



Substituting the general solution 

4>(x) = A cos kx + B sin kx 
back to the integral equation (3.36), we find the secular equation for the eigenvalues 



ka sin k(b — a) = cos k(b — a) 



with the solution 

g\N c +l) 



6tt 4 (2j + l) 2 iV c 



'ln5 N 
lni 



+ 



'ln5 N 
lni 



(3.41) 



(3.42) 



where j = 0, 1, 2, .... The corresponding eigenfunction reads 

<j>j{x) 



\ 



In 



- cos kj (x — b) 



which satisfies the orthonormal condition 

rb 



I dx<f>i(x)<j>j(x) = 5i 

J a 



(3.43) 



(3.44) 



The pairing temperature is determined by the condition that the largest eigenvalue, E , 
crosses one and we obtain that 



In 



kBTpair 



+ 5^2-^(1024^2^ 2 ) 



(3.45) 



which fixes the leading order term, the non-BCS exponent of (2.8), together with a part of 
the sub- leading contributions of the perturbation series (3.1). 

In appendix B, we shall derive the same eigenvalue condition (3.41) by summing up the 
perturbative expansion (3.10) of the Fredholm determinant for the equation (3.36). 
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D. Higher order corrections 



It follows from (3.42) and the condition E — 1 at T — T pair that 

ln-^-~- (3.46) 
k B T g 1 ; 

for the temperature around T pair . Therefore the order of the perturbation series (3.1) does 
not follows the explicit orders of the kernel diagrams. On the other hand, the forward 
singularity is not expected to be enhanced in higher order diagrams ( an example will be 
given in the next section ). The higher order kernel diagrams are not expected to contribute 
to the leading order term. In this subsection, we shall develop the systematics to handle the 
higher order corrections. 

To begin with, we construct a complete set of state vectors that diagonalizes the kernel 
(3.9) to the leading order with discrete Matsubara energies. We introduce 

<u,p\0> = C o f (v,p) (3.47) 
< ",p\j ^ > = Cjfj(u,p), 

and < u,p\a > with the adjoint expression 

(a48) 

where fo(v,p) is given by (3.35) with exactly O of (3.43) and fj{v,p) (j ^ 0) by the 
same ansatz (3.35) but with slightly rotated from <pj (3.43) to adapt to the ortho-normal 
condition with respect to discrete Matsubara energies. 

1 f°° 

<i \j>=15^1 dp<i\u,p><u,p\j >=Sij. (3.49) 

The complementary set of states, \a > is chosen such that < j\a >= and 

< a'\a >= 8 a i a . (3.50) 
which together with \j >'s make up a complete set of basis, i.e. 

^|j'Xj| + ^|a><tt| = l (3.51) 

3 a 

A clear cut zeroth order kernel operator is defined by 

K = Y,Ef\j>< 3 \ (3.52) 



with 

= ^: \ \:.^r (I"") - (3-53) 



.(Q) 9 2 (N c + l) ( l- 2 
6tt 4 (2j + lfN c \ e, 
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The contents of AK 


The contribution to b 


discrete v vs. continuous v 


7 + In 2 


correction to the approximation (3.33) 


none 


restoring exact au 


none 


restoring exact oe 


none 


cutoff 5 


none 


states below Dirac sea 


none 


higher order kernel diagrams 


-^(vr 2 +4)(iV c -l) 



TABLE I. The contributions to the sub- leading term of (3.1) 



We have then K \j >= Ej\j > and K \a >= 0. The difference between the exact kernel K 
and K Q , 

AK = K-K (3.54) 
generates the perturbation series of the eigenvalue, i.e. 

E l = Ef> + < j m\ j > + ^ <^ >< ^> (3.55) 



E? ! - E, 



v < j\AK\a >< a\AK\j > 



and the transition temperature is determined by E = 1. 

The analysis of perturbative corrections to the sub-leading order is rather laborious. We 
find few contributions from the first order term of (3.55) and none from the second order 
terms. Except for the contribution from the quark self-energy to be discussed in the next 
section, the details of the analysis will not be presented here. The interested reader is 
referred to the original work [16]. The result of the analysis is shown in Table I. 

A similar perturbation theory has been developed for the energy gap below T c [37]. 



E. Pairing at a nonzero angular momentum 

For equal helicity pairing, there is no net spin projection in the direction of the relative 
spatial momentum of the pairing quarks. The total angular momentum J > and the 
corresponding angular wave function is the Legendre polynomial, Pj(cos9). The partial 
wave component of the one-gluon exchange vertex is given by [16] 

7V,p'|M = -^(l + ^) j' i d(cos9)P J (cos9) (3.56) 

3 _ cos g_(l + cos g)| E |a! 1 + cosfl n 

X > " + \P- P? + mlf M {^) + \P ~ P? + mlf E (^) J ' 
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It follows from the forward singularity and the fact that Pj(l) = 1 that the leading order 
contribution should be independent of J. To figure J-dependence in the sub-leading order, 
we write 

P J (cos 9) = 1 + [Pj(cos 9) - 1] (3.57) 

and correspondingly 

Since the expression inside the bracket of (3.57) smears the forward singularity, the HDL self 
energies, <jm and <j£, can be dropped and the momenta p, p' can be set at \i when evaluating 
sj. We find that 

s J=n dcos6— — . 3.59 

2 J -i 1 — cos 9 

Using the recursion formula for Legendre polynomials, 

(J+ l)P J+ i(cos#) - (2J+ l)cos#Pj(cos#) + JPj_i(cos#) = 0, (3.60) 
we find that, 

(J + - (2 J + l)sj + Jsj^! = -<f J0 . (3.61) 
Since Co = 0, the solution to (3.61) for J > 1 is, 

*j = -E^ ( 3 - 62 ) 

n=l ' t 

For a cross helicity pairing, the spin projection in the direction of the relative momentum 
of the pairing quarks is one and the total angular momentum J > 1 (the orbital angular 
momentum is always perpendicular to the relative momentum). The corresponding angular 
wave function is the Wigner P-function df^Q) defined by 

d J M , M (9) =< JM'\e- lJ y \JM > (3.63) 

and the corresponding partial wave component of the vertex function reads [16] 

fV^IM = (l + ^) (2J+1) £ i d(cos9)d J 11 (9) (3.64) 



x 



l + cosfl-(l + cosfl)|-| ^ l_+cosfl 



\v - v>f + \p - p'P + m y M (^) \P-PV + mlf E (x) 
Following the same strategy for equal helicity pairing, we find that 
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FIG. 5. The quark self-energy diagram. 
d J n (9) cos 2 f 



[j = - dcostf^V \ 3.65 

J 2 J-i l-cos# v ; 

1/ J J + l \ 

= 2i Sj+ 2TTT Sj+1 + ^7TT Sj - 1 > 



The leading order degeneracy among different angular momentum makes pairing vulner- 
able to distortion in the presence of an anisotropic perturbation. This will impact on the 
LOFF pairing to be discussed in section V. 



IV. THE HIGHER ORDER KERNEL DIAGRAMS. 

A. The quark self-energy function 

Another important consequence of the forward singularity is the non Fermi liquid behav- 
ior of the quark self energy, which has been considered in solid state physics in the context 
of electrodynamics with quarks replaced by electrons and gluons by photons [36] . The kinks 
of the single particle occupation number distribution function at Fermi surface and T = 
is smeared and the electronic specific heat acquires a term proportional to T In T at low T. 
Both effects are too small to be observed in a nonrelativistic system. For an ultra relativistic 
system like the dense quark matter considered here, the impact is large and will lower the 
color superconductivity scales significantly because of the suppressed quasi-particle weight 
towards the Fermi surface. The quark self energy function S(P) is represented by Fig. 5, 
with the gluon line containing HDL resummation. Standard Feynman rules yield 

S(P) = C f Z(P) (4.1) 
with Cf — TfTj = ^j^ 1 for Tf in the fundamental representation of SU(N C ) and 

s(P) = -jT, j ^ys v ^ L M L + ( 4 - 2 ) 



where L = (I, —u), P = (p, —v) with uo and v the Masubara energy of the internal gluon 
and quark propagators. The bare quark propagator is given by 

S(P)=p (4-3) 

with P = 74 (fi + %v) — i'y • p. All gamma matrices are hermitian. Introduce the infrared 
sensitive region of the loop energy- momentum, < I < l c and — u c < uj < u c with l c , uj « /i, 
we have 
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S(P) = ~<(P) + ~>(P) 



(4.4) 



with the superfixes denoting the integration inside and outside the infrared sensitive region 
and we shall evaluate S < (P) only. 



~<(P) 



9 



2 r l 







4tt 2 



dll 2 / d(cos#)- J2 







\u>\<u>c 



£-i(v + v) 



8tt 3 



o 



dll dvV(l,v)F(u,p;l,v), 



where 



F(v,p;l,iv) 



rp-fi+i 

/ 



Carrying out the integral for p = fi, we obtain that 



P(i/, fi;l,u) — 2^74 tan 1 — \-2 U 



IV + V 



I 



7-p(l - 



and 



-F(^, m l, v) = 2m^6(u + {(j - l )2 - /2 



I 



74 



tan 



I 



iv + V 



2 

+ 7 



V + V -i 

2— —tan" 1 

v 



I 



+ 



2(u + v) 2 + l 2 



iv + v (iv + v) 2 + l 2 



1-p, 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



where the delta function comes from the discontinuity of the inverse tangent function. We 
find the energy dependence of the self-energy by differentiating, 



= g 2 [A(v) + 



(4.9) 



with 



and 



A(u) 



74 / dl 



An 2 Jo i2 + v 2 + m 2 D f M (=^ 



% rlc ru c , 

B{u) = — / dl duV(l,u){ 

2lt 6 JO J-L0 C 1 



il 



r74 



2— —tan" 1 

v 



u + z/) 2 + I 2 
I 2(v + u) 2 + I 2 



+ 



v + u (v + u) 2 + P 



1-p 



•} 



(4.10) 



(4.11) 



Noting that the asymptotic behavior, that /m(<z) ~ \ z \ for \z\ << 1, a scale Iq may be 
introduced to divide the integration in A(y) further into two: \v\ « Iq « (m 2 D \ v | ) . For 
I < Iq we have the contribution 
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Wo / 1 rk 1 1 Wo / 3 

/ dl y-r < —r / dl y-r < y-r / dll ~ « 1. (4.12) 

7o p + u 2 + m 2 DfM ^ mlJo / M (j) m 2 D f M (z)Jo m 2 D \v\ 

All of the inequalities follows straightforwardly from the definition of Iq except for the second, 
which is because of /m(^/0 being a monotonically decreasing function of /. Neglecting this 
sub-leading contribution, we find an infrared logarithm in A{y) arising from the second 
region (namely the region l < I < l c ). The integration B{y) is, however, finite in the limit 
v — > 0. We end up with 



d_„ <( , _j_ f l ° I 2 ig 2 4/ e 3 

Q V - ^ p >- 47r2j / /o dl p + v 2i + m l\ v \ ~ I2vr 274 nm 2 D \v\ 



and 



d 

— E < (i/,p) = ig 2 B(v) = finite (4.14) 



at p = fx. It follows then that 



^)U=-|^ /74 Wn-g^ + ... (4.15) 



B. The contribution of the quark self-energy to the sub-leading term 

Because of the trade off (3.46), the quark self-energy represents a order g correction to 
the leading order kernel diagrams and therefore contributes to the sub-leading term of (3.1). 
The calculation of this contribution proceeds as follows [14]: 

The quark propagator with the self-energy correction reads 

S(p, v) = % —— + AS(p, v) (4.16) 

iv — p + fj, 



where 



with 



A5(p, v) = -^E(P, v) (4-17) 



id 2 4/ 3 

E(p,v) = u(p)Z(P)u(p) = -T^ 2 Cfv\n—f-. (4.18) 
The corresponding correction to the kernel is 

2 

AK(v',p'\v,p) = ^ J (v',p'\v,p)[AS(p, v)S(p, -v) + S(p, v)AS(p, -v)\ (4.19) 

Applying the perturbation theory developed in the previous section, we find the first order 
shift of the maximum eigenvalue of the Fredholm kernel [14] 
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] _ _ /"CO POO _ 

^o = ^EE d P' dpf (u',p')AK(u',p'\u,p)Uu,p) 



(4.20) 



E (o) r dv_ roo dpp 2 f j {u p)[AS{p ^ v)S fa _ v) + 5(pj -i/)]/ J (i/,p) 

J-KkuT H JO 



"A) 



iV r 127T 2 J7rA: B T 7T 



5 r/z/ 4/ 3 



dp 



vrm^lz/l jo [z/ 2 — (p — ji) 



212 



e jo 



iV r 3tt 4 



<? 2 iv c 2 -i 

48tt 4 iV, 



7r 2 + 4)^ 0) In-, 



1 



where the wave function f^iy^p) is given by the ansatz (3.35) with O by (3.43). Upon 
setting the corrected maximum eigenvalue, 



En 



6vr 4 N c 
to one, we end up with 



(4.21) 



e 



7 - In 2 + 6sj + c ig (tt 2 + 4). 



iV c + l 9 
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(4.22) 



According to the arguments given in the rest of this section, we have completed the sub- 
leading term of (3.1) and the scaling formula of the pairing temperature (2.8) follows from 
(4.22). 



C. The infrared behavior of the vertex functions. 

The contribution of the quark self energy to the sub-leading term of the perturbation 
series (3.1) raises naturally the issues of the contribution from the vertex corrections, since 
they are intimately related through BRST identity of the nonAbelian gauge invariance. On 
the other hand, the Fermi sea makes the infrared behavior of the vertex function to depend 
subtly on the order of the infrared limit, i.e. the limit of zero energy transfer and zero mo- 
mentum transfer. Consequently, we find that the contribution from the vertex correction is 
beyond the sub-leading order without contradicting to the BRST identity. In this subsection 
we shall explore the infrared behavior of the vertex correction and demonstrate the absence 
of its contribution to the sub-leading terms of (3.1). The matching of the BRST identity to 
the leading infrared logarithm will be discussed in the next subsection. 

The vertex corrections are listed in Fig. 6. We begin with the first of them, which presents 
in an abelian gauge theory, say QED, as well. In order to simplify matters further, we shall 
put the spatial momenta of both external quarks on the Fermi surface, i.e. p — p' — fj,. We 
have 

Af)(P', P) = gTf^Tfk^P', P) = gT) {-^f- + C,) A M (P', P), (4.23) 



21 



where 

A,(P', P) = -/ll\ J2^ V ^ ^ S ( L + P ')l,S(L + P) 7 „ (4.24) 

and C a d5 c ' c = f abc ' f abc = N C 5 C ' C . This may be written in terms of two integrals, one inside 
and one outside the infra-red sensitive region: < I < l c , —u c < u < uj c with l c , u c <C fi, 

A M (P', p) = [a (o)< (P', P) + A (a)> (P', P)] . (4.25) 

with P = (— ip, 1). As we are only interested in the leading infra-red behavior, we evaluate 
A(a)<(P' ; P) on iy ; 

A (o)< (p/ ) p n _ 9 f c dU 2 f d ^ cos ^ , _^.p cos 2 ff\\(p^ p>- L) sin 2 9, (4.26) 
8ir d Jo J -i 

— # 

where ( = v + i£, £ = \l + p \ — fi and (' and £' refer to v\ p' . The logarithm in (4.27) 
introduces a branch cut which we may take to lie along the positive real axis and the contour 
to run above and below it in the normal fashion. As shown for the self-energy, it is only the 
discontinuities that occur as poles cut the contour and branch cut that induce the infra-red 
singularity. Hence we need only focus upon the second and fourth terms in (4.27), since 
the other terms are regular. Using the convention that arg(— oo) = along the negative real 
axis, we find that the contribution of these poles reads 

A disc .(P', P; L) = -L- [Sign(OX>(J, C) - Sign(£')^, C')] , (4-28) 

where the sign function comes from the discontinuity of ln(— u) crossing the cut. 

We are now in a position to examine the limit P^ — > P^ and we shall consider the two 
different orderings of the limits in turn: 

i) lim \imA disc .(P', P; L) = -7rSign(o|-P(Z, C), (4.29) 
ii) lim. lim A disc .(P', P; L) = in-?- [Siga^)V(l, ()] . (4.30) 

p'^pv'^v (jq 

In both cases we are looking at the infra-red limit, and thus fix the external momentum to 
be p = p' = fi. 

First of all, considering case i), changing the integration variable from dcosO to d£, it is 
straightforward to find, 

lim Km A«< (P',P) = ~f~ 2 f° dll f Sign(o|-Z>(', C) (74 - ^7 • p) (4.31) 

u'~-*up'^p p=li ^ti z Jo J— i av v L z ' 

2 r l 



74 l C dll[V{l,u) -V(l,v + il)] 
Jo 



W 
4tt 2 

ia 2 4/ 3 
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where the second term inside the brackets of the second line of eq.(4.32) is free from loga- 
rithmic divergence in the limit v — > and is denoted by ellipses in (4.32). 

Secondly, considering case ii), we find that differentiation gives two terms which will 
cancel in the leading order, 

Sign(O^Z>(J, + 26(t)V(l,Q 

(4.33) 

The first term is identical to that evaluated for case i) above. With the same approximation, 
the second term is, 

zl r lc d i 12 ~ zi i n l i 

An 2 Jo P + m 2 D u Yin 2 m 2 D v 

The leading contributions of the two terms cancel and in this ordering of limits the vertex 
is free from the infrared logarithm. 

When this vertex diagram is inserted in to the kernel of the Dyson-Schwinger equation 
for di-quark scattering, the important contribution to the pairing, as is determined by the 
forward singularity of the one gluon exchange, comes from 

\p -pf ~ m 2 D \v' - v \ ( 4 - 35 ) 

which means that 

\p — p\ » \v' — v 

according to (3.27) and is away from the region of energy-momentum transfer, 

\p-p\ « W -v\ (4.37) 

in which the vertex correction is enhanced logarithmically. This is the reason for the ab- 
sence the absence of the vertex correction from Fig. 6a in the sub-leading terms of (3.1) and 
the statement was also verified numerically [17]. We have also checked explicitly that the 
vertex correction Fig. 6b does not contribute to the sub-leading term. The absence of the 
contribution of Fig. 6c to the sub-leading term follows from the BRST identity discussed 
below. 



lirn Jim A^<(P', P) = |4 t dl I f d^U - if 
j'-»p 8tt 2 Jo J -i v r 



1-P 



(4.34) 



(4.36) 



D. Matching the BRST identity. 

In this subsection, we shall show that the absence of the vertex corrections in the sub- 
leading terms is not in conflict with the BRST identity of the gauge symmetry. The BRST 
identity is the generalization of the Ward identity of an Abelian gauge theory to a non- 
Abelian gauge theory and it takes the form for one-loop diagrams incorporating HDL gluon 
propagators [17] 

(P' - Pyh^P', P) = Tj(£(P') - S(P)) + (P' - P)^(P', P). (4.38) 
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(a) (b) (c) 



FIG. 6. The physical radiative corrections to the quark-gluon vertex; a) A p , the abelian 
vertex, b) A 1 ^ , the tri-gluon vertex and c) Aj^ , the triangular vertex. 




(a) (b) (c) 

FIG. 7. The non-physical ghost diagrams generated by the BRST transformations; representing 
a) R l ( a \ b) and c) R l{ - C \ 



The physical quark-gluon vertices = A 1 ^ + A 1 ^ + A 1 ^ are represented in Fig 6. The non- 
physical ghost-quark vertices generated by the BRST transformation, R l = R 1 ^ + R 1 ^ + 
R l<yC \ are represented in Fig. 7 and vanish for on-shell Minkowski momenta P and P' at 
H = 0. 

To prove this identity, we start with the Feynman amplitudes of the quark self energy 
diagram of Fig 5. 

a 2 r d 3 l 

S(P) = -C f 9 -Y^J —V^L^P + L) lu . (4.39) 
and the vertex diagram of Fig 6a. 

^ a \P', P) = -TfT\Tf^ £ / ^V vp (l \u n \) lv S(P' + Lh,S(P + L) lp . (4.40) 

Using the standard trick, 

(P' - PyS{P' + L) ltx S(P + L) = S(P' + L) - S(P + L), (4.41) 
we may trivially rewrite this expression in terms of the self-energy, 

(P' ~ PYAf\P',P) = T l f (l- ^ [E(P') - E(P)]. (4.42) 

For QED C ad = and (4.42) becomes the ordinary Ward identity. However, for non-abelian 
gauge theories, the group theoretic factor do not match but we have additional vertices 
which cancel the extra term. 

The Feynman amplitude of the second vertex diagram of Fig. 6 reads 
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A W(P',P) =frnn T rn T nl^ J ^L Vl , X {L){-i)T » Xp {L, L - Q)V pl/ ,{L) lv S{P + L) 7 „, 

(4.43) 

where if lmn T™Tj = — and r MAp (P',P) denotes the bare tri-gluon vertex. Using 

the explicit expression of Y il x p {P l , P) and the Dyson- Schwinger equation for a HDL gluon 
propagator, 

V^{K) = D, V {K) - iV, p {K)Tl pX {K)D Xv {K) (4.44) 
with D(K) the free gluon propagator. We find that 

(P' - Pyv v , x {P')T» Xp {P\ P)V pu {P) = ig[V^(P', P) + V^(P', P) + V^(P', P)], (4.45) 
where 

V^(P', P) = i \V V , V {P) - V V , V {P')\ , (4.46) 
KS(P', P) = V V , X {P>) [U Xp (P') - U Xp (P)} V pv {P\ 
V^(P', P) = A(P')P^P' x V Xu (P) - V V , X {P')P X P V /±{P). 

with A(P) = — ^ the ghost propagator. Correspondingly, we have 

(P' - Pykf \P', P) = I[(P', P) + I l 2 (P', P) + I l 3 (P', P) (4.47) 

with 

fyP', P) = -CadT) J ^V$(L, L-P' + P) lv S{P + L) 7 „. (4.48) 

It is straightforward to show that, 

(P'-P)"A^(P',P) + /{(P',P) = #Tj[£(P') - £(P)] (4.49) 

and I l 3 (P',P) = (P' - P) M P| t (P',P). The second term on RHS of (4.47) is to be canceled 
by the third vertex diagrams, Fig. 6c. 



Denoting by r'™"(P', P) the triangular vertex in Fig. 8 and using the identity 



we find that 



P' ~ P)^X( P> > P ) = wf lmn in Xp (P') - U Xp (P)} (4.50) 



(P' - Py\ l jf\P', P) + I l 2 (P', P) = (4.51) 

and the BRST identity (4.38) is established. 

As is shown in the previous subsection, the infrared logarithm shows up in the limit 
v 1 — > v following j? — > p, i.e. 



25 



FIG. 8. Quark loop with three external gluons, T 1 ^. 



lim lim kf\P', P) = T\ ( -Q* + C f ) In — ^-"T, (4.52) 

In the same order of limit, the BRST identity becomes 

lim limjAi(P', P) - R[(P f , P)] = - gT l ®Z(P). (4.53) 

It has been shown in [17] that both Af\P',P) and R l 4 (P',P) are free from the infrared 
logarithm in this limit, but A l (°\P', P) contributes one, 

lim Inn A l j c \p', P) = i^—C adl ff) In ~^ n , (4.54) 
v'^up'^p 247T 2 ' 7cm z D \v\ 

which together with (4.52) and (4.53) match the leading logarithm of the BRST identity. 

The presence of Aj c )(P', P) in the BRST identity is not surprising in terms of the stan- 
dard proof of Ward identity by summing up all possible insertions of a gauge boson line into 
a fermion self energy diagram. Since HDL resummation is employed for gluon propagators, 
the inserted gluon line may land on one of HDL's and give rise to A 1 ^. 

It worth mentioning that the nonzero chemical potential makes the triangle vertex of 
Fig. 8 non-vanishing for QED as well. But it will not contribute to the Ward identity with 
HDL photons because the RHS of (4.50) vanishes. 

The BRST idensity in the superconducting phase is dicussed in [20]. 



E. The crossed box diagram with two gluon exchange. 



Now we come to the last 2PI vertex shown in Fig. 4, which contains two HDL gluon 
lines crossing each other. This was estimated in [8] by a renormalization group argument. 
The explicit calculation of this diagram in [14] will be outlined below. 

Apart from a group theoretic factor, the contribution of the crossed box diagram to 
7+ +++ reads 

xW)l r S F (P + T,^u(v)\H-v')l,'S F (-P + (>) vll (-p)], (4.55) 

where \p\ = \p'\ — /i, P — |(p +p'), q = p — p' and the discrete Matsubara energies have 
been replaced by continuous Euclidean energies with the external ones set to zero. Ignoring 
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the Coulomb propagator, the contribution from the small scattering angle, \6\ < 9 <C 1, 
and the infrared region |/| << p, \u\ « /i, 

n 1 4 f 9 ° r/i n f d0J d 3 l -> q . ,r* Q . 

= - 2 9 J_ eo dds ™ d J m to J2^ V ^ 1 - 2^ V ^ 1 + 2'^ 

x[u(P )liSF{P + l,iu)^ju(p)][u(-p )^ r S F (-P + l,iu)ji>u(-p)], (4.56) 
is bounded: 5j R < /, where 

J =327vi) M L df>dh \r% + k\ P \){A+ L|p|)(p 2 + £|)(p 2 + El)' (4 ' 5?) 

Here E± = \P ± Z — 1, r± = |/ ± (f|/p, p = M/p and /t = f-^r- Transforming the 
integration variables from 9, r to E±, r±, we end up with / = 32 ^ 4 ^ 2 Jo dpK(p) where 

K( P ) = J dE + dE_drldrl J ^ - 4^^^^ + (4-58) 
with the Jacobian 

J = [(E+ - EJf - 4(r 2 + r 2 _)(E + - EJf - 16(£+ + £_) 2 + 16r 2 r 2 ]" 1/2 . (4.59) 

As p — > 0, we find that K(p) — > const xp~ 2 / 3 and / = finite. The absence of the forward 
singularity for the two gluon scattering vertex renders its contribution to the transition tem- 
perature smaller than the leading term by a factor of g 2 / In -^f - Therefore its contribution 
is beyond the sub-sub- leading order according to the trade off (3.46). 



V. LOFF PAIRING WITH A MISMATCHED FERMI SEA OF QUARKS. 



A. A mismatched Fermi sea of quarks. 



In the previous two sections, we took the massless limit of quarks, which is a good 
approximation at ultra-high baryon density. As the density is lowered, the nonzero masses 
of quarks have to be taken into account and different flavors will not have the same Fermi 
momentum. The phase space available for pairing will be reduced and new orders may arise. 
Crystalline color superconductivity discussed in [38] is a potential candidate. 

In order to see how the Fermi sea of quarks is mismatched because of nonzero masses, 
consider the equilibrium of u, d, s quarks and electrons. The number density of each flavor 
is 

A; 3 

n, = ^ (5.1) 

where kf the Fermi momentum of each flavor with / = u,d,s. The number density of 
electrons is n e = with k e the Fermi momentum of electrons. The extra factor 3 for rif 
comes from the three colors. Ignoring the interactions, the total energy at T = is 
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3 r^ u i 3 r^d i 3 f^ s i 

E = — dpp 2 \l p 2 + ml H — -/ dpp 2 J p 2 + m\ H — -/ dpp 2 Jp 2 + m 2 (5.2) 

TT Z JO v 7T JO 7T Z JO V 



3 /■*« , , ~ 3 r fc d 

^io 

The total baryon number of the system is 



b = ^(n u + n d + n s ) (5.3) 

and the electric neutrality reads 

2 11 . . 

-n u - -n d - -n s - n e = 0. (5.4) 

The Fermi momenta kf and fc e are determined by minimizing the energy (5.2) subject to 
the constraints (5.3) and (5.4), 

dF dF . . 

where 

1 2 11 

F = E - -fi B {n u + n d + n s ) - /i Q (-n u - -ra d - -n s - n e ) (5.6) 

with n B and /iq two Lagrangian multiplier. Under the approximation kf » rrif, we end 
up with slightly different Fermi momenta for different flavors. 

k u = k F (5.7) 



2 ™2 



k d = k F + 



mz — m 



d 



Ak F 

2 

777 . - 

k s = k F + 



4A; F 

2 1 

where /cp = 7T3&3 is the universal Fermi momentum in the chiral limit. As the baryon 
density is lowered further such that k F ~ m s , more electrons will be involved for the electric 
neutrality. 



B. Pairing instability at a nonzero total momentum. 

A mismatched Fermi sea is known in a metallic superconductor with ferromagnetic impu- 
rities. The exchange interaction between the electron spins and the impurity spins displaces 
the Fermi momentum of each spin while levels up the total energy (kinetic and exchange) 
of each electron on the Fermi surface. The phase space available for BCS pairing (pairing 
between the electrons of momentum ±p) is reduced and other pairing states may arise. A 
potential candidate of them is that suggested by Larkins, Ovchinnikov, Fudde and Ferrell 
(LOFF) which pairs the electrons of momenta not exactly equal and opposite to each other 
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FIG. 9. The phase diagram showing the competition between BCS pairing and LOFF pairing. 

[39]. The resultant Cooper pair in this case carries a net momentum, 2q, and the gap 
parameter supports a crystalline structure in the coordinate space, 

A(f) = A e 2 ^. (5.8) 

Denoting the difference between the Fermi momenta of the two pairing quarks by 26, the 
phase diagram on the T — S plane that shows the competition between BCS pairing and 
LOFF pairing is displayed in Fig. 9. At zero temperature, BCS state persists for < S < S 
and LOFF state wins for S < S < S nmx . Along the line of superconducting transition, it 
is from the normal phase to BCS phase for < S < Si and it is from the normal phase to 
LOFF phase for Si < S < S max . No long range orders exist for S beyond 5 max - The BCS 
phase and the LOFF phase below the transition temperature are separated by a first order 
phase transition. 

The Dyson-Schwinger equation of section II can be easily generalized to the diquark 
scattering at a nonzero total momentum, i.e. 

r<r, s (Pf\Pi) = r q -AP f \Pi) + ^ E / -^K^{p f \p)r 9tS (p\p^ (5.9) 

where the kernel is given by 

K^{P'\P) = r<rAP'\P)S(Q + P^ + S)S{Q -P,n-8) (5.10) 

with f the 2PI part of the scattering vertex. The letter P's in (5.9) stand for the relative 
four momenta between the pairing quarks with P = (p, —v) and the letter Q for their 
average four momentum with Q = (q, 0) (total momentum being 2Q). The color-flavor and 
spinor indices have been suppressed, and the Fermi momentum of each quark propagator is 
indicated explicitly. The eigenvalue problem of the kernel reads 
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The maximum eigenvalue is a function of T, fi, 5 and g and the pairing instability is given 
by the condition 

E max .(T, f ji,6,q) = l, (5.12) 

which gives the pairing temperature as a function of /i, 5 and q, i.e. T pair (/i, 5, q). The 
transition temperature for a given /x and 5 is determined by 

T c = max(T pair (/i,<5,g)|Vg). (5.13) 

If this occurs at q — 0, the superconducting state below T c is of BCS type. If this occurs at 
q 7^ the superconducting state below T c is of LOFF type. 

In the next subsection, we shall tackle the eigenvalue problem (5.11) for a point inter- 
action and reproduce the results known to condensed matter physicist [40]. The solution 
for QCD will be presented in the last subsection, where we shall show how the forward 
singularity of the one-gluon exchange widens the LOFF window in the phase diagram on 
T — 5 plane [41] [42] and distorts the shape of the order parameter below the transition 
[42]. 

C. LOFF pairing for a point interaction. 

Consider the Hamiltonian with two species of fermions with different Fermi momenta 
interacting with each other via a point interaction, i.e. 

H = YNp - ^ - 6 ) a hp + (p-v + - ^ £ GfWq+^l-fk-^g+P ( 5 -i 4 ) 

p f,p,q 
with Gp£ = G9(uo — \p' — h\)6{ujq — \p — and G > 0. The 2PI vertex in this case is 

t^{P'\P) = -G6(uj - \p' - (i\)6(uj -\p-fjL\) (5.15) 

and the kernel (5.10) reads 

K^(P'\P) = -G9(lu - \p' - pl\)9(u; - \p - pi\)S(Q + P, fi + 5)S(Q -P,n~5) (5.16) 

with 

S(P,im) = - (5.17) 

iv — p + jj, 

Only the zero angular momentum channel is present for pairing. 

The integral equation (5.11) is completely separable and admits an exact solution. There 
is only one nonzero eigenvalue given by 
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E = 1 + 



111 A () 



In 



+7 + V ; ( _ ) ~~ — Imln ^ 2 



i(S-g) N 

2 2nk B T / 



(5.18) 



where A is the energy gap at T = 0, 5 = and g = 0. The corresponding eigenfunction is 



f(P) 



1 



4.>tMr^-iP-.i). 



(5.19) 



The eigenvalue problem (5.11) can also be solved by the perturbation theory developed in 
section III with the kernel of the zeroth order given by (5.16) at 5 = and q = and the 
perturbation by 

A K{P'\P) = K S ^P'\P) - K ofi (P'\P) (5.20) 
= -G6(u - \p' - n\)e(u - |p - fi\)[S(Q + P,fi + 5)S(Q -P^-8)- S(P, n)S(-P, //)]. 

Eq.(5.19) is the eigenfunction of both the zeroth order kernel and the perturbation. Thus 
the first order perturbation gives rise to the exact eigenvalue (5.18). 

The pairing instability is located according to the equation E max = 1, i.e. 



In 



Ap 

nk B T 



+ 7 + ~~ — Imln 



1 _ i(S-q) 

2 2wk B T 



v (l _ i(S+q) \ 
L \2 2wk B T J 



o, 



(5.21) 



with ip(z) = dln ^} z " > and the transition temperature is determined by (5.13) [40]. The 
numerical solution to (5.21) yields the values of relevant parameters in the phase diagram 
Fig. 9. 



tfi = 0.606A 
<Wc = 0.745A . 



(5.22) 



The BCS state amounts to pair quarks of equal and opposite momenta at disadvantage 
of the mismatched Fermi sea. The quasi particle operator is related to ap and bp by a 
Bogoliubov transformation 



Oit. 



Up Vp 
Vp Up 



ap 



with 



1 r 



U- 



1 + 



p-fl 



1 r 



(5.23) 



(5.24) 



p - 11) 2 + A 2 

The ground state, |BCS >, is annihilated by a or f3 of an arbitrary momentum. The quasi 
particle energy 
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E} = ±5 + y/ip - ^ + A 2 (5.25) 

opens a gap A — 5 for 5 < A. 

The LOFF state amounts to pair the quark of momentum q + p with that of q — p and 
the quasi particle operator is given by a modified Bogoliubov transformation: 

(5.26) 

q-P / 

The quasi particle energy is 

Ep = ±(q cos 9 + 5) + y/(p - /i) 2 + A 2 (5.27) 

with # the angle between p and <f and can be negative. The corresponding ground state 
reads 

iloff >= n 4 n 4\ > ( 5 - 28 ) 

E+<0 E p <0 

with the state | > annihilated by a and (3 for all p. Upon comparison between the free 
energy of a BCS state and that of a LOFF state at T = 0, we find that BCS state extends to 
5 ~ ^ and is replaced by LOFF states up to 5 max - The BCS state and the LOFF state are 
separated by a first-order phase transition. Using Ginzburg-Landau approximation to the 
gap equation, the authors of [43] argued that the transition from the normal phase to LOFF 
phase is also of the first order because of the negative quartic term for certain crystalline 
structures of the LOFF state. 



D. LOFF pairing with one gluon exchange. 

Consider an idealized problem of two flavor pairing, Nf — 2, with a Fermi- momentum 
difference 25. The vertex function of one-gluon exchange within the color anti-symmetric 
channel and with a total spatial momentum 2q, T^s{P'\P), contains all angular momentum 
channels and reduces to T A of (3.7) at q = and 5 = 0. The corresponding kernel of 
Dyson-Schwinger equation for di-quark scattering, eq. (5.10), is diagonalized perturbatively 
in sections III and IV under the same condition. The maximum eigenvalue within each 
angular momentum channel given by 

EL,. = 1 + ilrC ^ ^ + 7 + Qsj) + 0(ln" 2 \), (5.29) 

where 

g 5 k B T c 



2567r 3 /x' 



(5.30) 



with T c referring to the transition temperature at 5 = 0, A is the energy gap at T = 
and 5 = 0, and E^ ax . here is the same as E of eq.(4.21) for Nf = 2. The corresponding 
eigenfunction to the leading order reads 
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/jAf(P) = Caiii( Ane7 ) (5.31) 

with Yjm(8,<P) the ordinary spherical harmonics. The adjoint of (5.31) is given by 

fMP) = (5.32, 
and the constant C of (5.31) is determined by the normalization condition 



I?/ 



!/ /? = 1- (5-33) 



(2^) 3 

Near the transition, In ^ = 0(g~ l ). We notice that the dependence on the angular momen- 
tum is of the sub-leading order. In another word, maximum eigenvalues of different angular 
momenta are all degenerate to the leading order as a consequence of the forward singularity. 

We shall treat the case with q ^ and 5^0 with the aid of the perturbation method 
developed in sections III and IV. The perturbing kernel is 

AK(P'\P) = K q - s (P'\P) - K 0:0 (P'\P) (5.34) 
= [f^(P'\P) -r 0fl (P'\P)]S(P,ii)S(-P,ii) 
= f ,o(P'\P) [S(Q + P,l* + S)S(Q -P,im-S)- S(P, n)S(-P, v)] 
= [f^(P'|P) - fo,o(P'|P)] [S(Q + P, ii + 6)S(Q -P^-5)- S(P, fi)S(-P, //)]. 

Since the relevant momentum scale of the vertex function is fi, the difference [T^ t s(P'\P) — 
f fi(P'\P)] is expected to be higher than r 0:0 (P'\P) by the order of 0(J^). The only 
contribution to the sub-leading term comes from the propagator difference [S(Q + P, \i + 
S)S(Q — P,/j,) — S(P,n — 5)S(—P,fj,)], similar to the case of a point interaction. If q = 0, 
rotation symmetry is maintained. The shift of the eigenvalue is simply the expectation value 
of the perturbing kernel, like the case of the point interaction and we end up with 



6E 



In! 



4)-4-*db)l' (5 - 35) 



where the second term inside the bracket corresponds to the limit q — > of the logarithm 
of the ratio of gamma functions in (5.18). If q ^ 0, the dependence of the perturbing 
kernel on the angle between the relative momentum and the total momentum breaks the 
rotational symmetry. Different angular momentum channels will be mixed and a degenerate 
perturbation theory has to be employed to figure out SE. 

In general, the maximum eigenvalue of the kernel (5.16) can be written as 

2 A 

E max = 1 + ( In ° + 7 + Pmax.) + sub-sub-leading terms (5.36) 
ln iV nk B T ' 

e D 

with p max . the maximum eigenvalue of the operator 

^-op. = s op. + fop.- (5.37) 
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One-gluon exchange 


Point interaction 


Upper limit of LOFF window 


<W = 0.968A 


<W = 0.754A 


Onset of LOFF along T c 


<5i = 0.606A 


<5i = O.6O6A0 


Oneset of LOFF along T = 


So < 0.707A 


S' Q = 0.707A 


Shape of the gap 


non-spherical 


spherical 



TABLE II. The parameters underlying a LOFF pairing. 



The operator s op , is diagonal in the angular momentum representation, 

< J', M' I s^. I JM >= 6sj8j,j8m>m, (5.38) 
while the operator v op , is diagonal in the coordinate representation (angle representation), 



v, 



W g/ S / ^/K^P-IPJ/^P) (5.39) 
, /1\ rl S — qcosO 



op. 



In the coordinate representation, the eigenvalue equation h op u = pu takes the form [42] 



dx 



,U{X)-U{x') f ( /l\ „ , fl .^-9^11 / \ / \ / r A r\\ 

{^(9) - Re ^k "'^FtJH x) = p u ( x ) ( 5 - 4 °) 



\x — x'\ 



2 2nk B T 



with x = cos^. For q 7^ 0, this integral equation can only be solved numerically and we 
obtain the values of the parameters in the phase diagram Fig. 9, i.e. 

5i = 0.606A (5.41) 
<W = 0.968A . 

Our result for 5 max is different from that of Ref. [41], which claims that 5 max = 00. The 
eigenfunction -u(cos#) corresponding to <5 max . is plotted in Fig. 10. This is to replace the 
spherical harmonics in (5.31) for M = for the zeroth order eigenfunction of the degenerate 
perturbation theory. 

The widening of the LOFF window comes from the angular momentum mixing, which is 
another consequence of the forward singularity. Furthermore the gap parameter right below 
the transition line follows the shape of the eigenfunction u(x) and is no longer spherical for 
5 > Si as is shown in Fig. 10. Because of the flexibility of the shape of the gap parameter 
for LOFF pairing via the angular momentum mixing, the other end of the LOFF window 
at T = is expected to be lower than that for a point interaction, i.e. S < 

The parameters underlying the LOFF pairing for one-gluon exchange and that for a 
point interaction are summarized in Table II for comparison. 
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FIG. 10. The eigenfunction u(x) corresponding to the largest eigenvalue of Eq. (5.40). 

VI. GINZBURG LANDAU THEORY OF COLOR SUPERCONDUCTIVITY. 

A. The general structure. 

Like a metallic superconductor, the physics of a color superconductor near the critical 
temperature can be described in terms of Ginzburg-Landau free energy [44]. Ginzburg- 
Landau theory is also a powerful tool to explore an inhomogeneous condensate as induced 
by an external field or by a nontrivial boundary condition, e. g. a vortex filament or a domain 
wall. Though derivable from the first-principle QCD Lagrangian in the weak coupling limit, 
the validity of the Ginzburg-Landau theory is not restricted within the perturbative region 
if the coefficients involved are regarded phenomenological 

The symmetry of QCD with three flavors u, d and s in the chiral limit is SU(3) C x 
SU(3)r x SU(3)l x U(1)b with the ordinary electromagnetic U(l) a subgroup of the flavor 
SU(3) i.e. C/(l) em C SU(3)r x SU(3)l- As was analyzed in previous sections, the most 
favorable pairing configuration is between two quarks of equal helicity at zero angular mo- 
mentum. The corresponding order parameter ^r(^l) for right (left) handed Cooper pairs 
carries the color-flavor indices of each pairing quark, i.e R(L)) C f\% an d is symmetric under 
a simultaneous interchange of them, i. e. 

(**L))J8 = (**l))SR- (e- 1 ) 

Since the pairing channel is 3, the order parameter is expected to be approximately anti- 
symmetric with respect to interchanging the two color indices, i.e. 

(*B(L))%% * -^R(L))%%- (6-2) 
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It follows from (6.1) that 

* -(*RW)%% (6-3) 

A small color-symmetric (sextet) component of the order parameter can be induced but its 
contribution to the free energy is of higher order as is analyzed in the appendix C. Dropping 
the color symmetric component, the most general form of the Ginzburg-Landau free energy 
functional consistent with the symmetry is 

r = J <?r^F\ 3 F\ 3 + 1(V x Af + ±Tr[(D* R )t(D9 R ) + (D* L )t (D* L )] (6.4) 
+ ^aTr[^ R + tft^j + I &T r[(^ R ) 2 + (^ L ) 2 ] 
+ ±b'[(Tr* R y R ) 2 + (Tr^ L ) 2 ] + C Tr(*^)Tr(^ L )]} 
where the gauge covariant derivative of the diquark condensate \1> reads 

= ~ igfr* ~ w^ c '^ R{L) )t h (6.5) 

-ie(q fl +q h )A(V R{L) ) c f \% 

with A = A l T l the color gauge potential, A the electromagnetic one and the trace is defined 
by TrMiV = M^^NJ^. Parity violating terms are neglected in (6.4) and an even parity 
condensate, 

V R = ^ L EE 

will be considered in subsequent discussions. The variational minimum of (6.4) with respect 
to the gauge potentials and the di-quark condensates satisfies a set nonlinear Ginzburg- 
Landau equations whose solution gives rise to the equilibrium configuration of the mean 
field \1>, A and A. 
On writing 

n% = ^^hhf^f (6-6) 

with $ a 3 x 3 complex matrix supporting the representation 1 © 8 under a simultaneous 
color-flavor rotation, the Ginzburg-Landau free energy becomes 

r = [ d 3 f\^F^F^ + -(VxI) 2 + 4tr (D<5>y (D$) + 4atr$ f $ + 6 1 tr($ t $) 2 + 6 2 (tr($ t $)) 2 ] , 

(6.7) 

where h = b, b 2 = b + 86' + 8c, 

/3$ = V$ - igrff 1 ® - iqAQQ (6.8) 
with T l the generator in 3 and the charge matrix 

Q = diag(-,--,--) = --^. (6.9) 
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The trace, tr(...), refers now to 3 x 3 matrices and we have rearranged the 577(3) generator 
to make T 8 proportional to the charge matrix. 

In the weak coupling limit, the Ginzberg-Landau coefficients can be derived from QCD 
one-gluon exchange or NJL effective Lagrangian and we have [45] 



a = ^ k BT c (T-T c ) (6.10) 
b' = c = 0. 

For a homogeneous color superconductor, eq. (6.7) becomes 

T = ft^atr^ + 6itr($ t $) 2 + 6 2 (tr($ t $)) 2 ]. (6.11) 

Depending on the coefficients of the quartic term, there are three different cases of the 
variational minimum of (6.11) [45] [29]: 

1) bi > and b\ + 36 2 > 0: The minimization of (6.11) produces 

$ = 0cfl£/ (6.12) 

with 



*» - t (6 - 13) 

and U an arbitrary unitary matrix. The order parameter supports the unit representation 
of the simultaneous color- flavor rotation and is therefore color-flavor locked. For U = 1, we 
obtain the standard form of the order parameter (2.22). with 4> s = and 4>a = 0cfl- The 
minimum free energy density reads 

r min . 12a 2 12/x 2 ^ ( Tc -T ^2 

where the last equality follows from the substitution of the weak coupling coefficients (6.10). 
We shall consider this region of parameters only for the rest of the section. 

Although a homogeneous CFL condensate can always be transformed to the standard 
form with U = 1 by a symmetry operation of (2.19), it is no longer the case in the presence 
of a vortex filament, which implement a nontrivial mapping of tti(U) [29]. 
2) ^ < and &i + b 2 > 

The variational solution in this case reads 

$ = diag(0, 0, Is e ia ) (6.15) 

with 

* s - fS? (6 ' ie) 
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and the minimum free energy is 



Tnim. 4a . . 

— = -b^r, (fU7) 

This case corresponds to an isoscalar [45]. Notice that the trace of (6.15) may be transformed 
away upon a flavor SU (3) rotation, e.g. $ — > $[/ with 



U= -1 . (6.18) 




Therefore this condensate belongs to the octet under a simultaneous color-flavor rotation. 
3) Other regions of the quartic coefficients: 

The stationary point of T becomes a saddle points and higher order terms of Ginzburg- 
Landau free energy have to be restored to stabilize the condensate. 



B. Characteristic lengths of a color superconductor. 

To identify various length scale that characterize an inhomogeneous condensate, we con- 
sider small fluctuations about the homogeneous one (6.11). We parametrize the order pa- 
rameter by 

* = fe L + ^(^) + f.(^) (6,9) 

with X's and Y's real, and form linear combinations of the ordinary electromagnetic gauge 
potential with the eighth component of the colour gauge potential 

V = A 8 cos 9 + Asm 9 (6.20) 

V = -A s sin9 + ^cos6», 

where the 'Weinberg angle' is given by 

tan0 = --^-, (6.21) 

Next, the free energy (6.7) is expanded to quadratic order in X's F's A's and A. We find 
that 

T = r min + J d 3 r{±(V x V) 2 + tr[(V x W) 2 + m 2 w W 2 ] lj[(Vx Z) 2 + m\Z 2 ] (6.22) 
+ \[{VX) 2 + m 2 H X 2 ] + l -[{VX l ){VX l )+m' 2 X l X l ] + i(W) 2 }, 

where W = ELi ? l W l and 
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\yi = A l -—VY l for Z = 1,2, - • -7 (6.23) 
m w 

Z = V - — VF 8 , 

m z 

The masses of the excitations provide us with the relevant length scales which are the 
coherence lengths, £ and £' defined by: 

2 2 
m | = (6 X + 3b 2 )(p 2 CFL = — , m'lj = &i0cfl = p , ( 6 - 24 ) 

that indicate the distances over which the diquark condensate varies and the magnetic 
penetration depths, S and 5' by 

m z = Ag 2 (p 2 CFL sec 2 e = , = m|cos 2 9 = . (6.25) 

The Ginzburg-Landau parameter that classifies the types of the response of the supercon- 
ductor to an external magnetic field is defined by the ratio 

k = | (6.26) 

and we shall see in the next subsection that the critical value that separates the type I 
superconductivity and the type II one is different from that of an metallic superconductor 
because of the multiplicity of the order parameter. 

The fluctuations considered so far are all of even parity, i.e. the relation = is 
maintained. The excitations that violate this relation are of odd parity and correspond to 
the Goldstone bosons associated with the chiral symmetry breaking induced by CFL. 



C. The magnetic response and the types of the color superconductivity. 

Before discussing the inhomogeneous condensation of a domain wall, we shall first de- 
termine the thermodynamical critical magnetic field of a homogeneous condensation with 
color-flavor locking. 

Like an metallic superconductor, the thermal equilibrium in a constant external magnetic 

— * 

field H is determined by Gibbs free energy, which is a Legendre transformation of (6.7), i.e. 

f = T - H ■ J d 3 rV x ^=fi + f 2 , (6.27) 

where 

fj = f d 3 f\j F ij F ij + x V) 2 + 4tr(z3$) t (z3$) + 4atr$ f $ (6.28) 
J 4 1=1 2 

+ & 1 tr($ t $) 2 + 6 2 (tr($ t <l>)) 2 -iJ- V x V sin6 

and 
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Super phase 


Normal phase 


$ 

V xV 

fl 


<ACFL 




p 12a 2 
S 4 6i+36 2 




Fsin# 


-Sl±# 2 sin 2 fl 



TABLE III. The comparison between the super phase and the normal phase in an external 
magnetic held. 



err 



\$ x V) 2 



H ■ V x Vcos# 



(6.29) 



The minimization of T 2 with respect V yields V x V = HcosQ and T 2 = — \VtH 2 cos 2 6. 
Since V corresponds to the unbroken U(l) gauge symmetry and does not couple to the 
order parameter, we shall concentrate our attention to Ti from now on. The results of 
minimization of T\ in the super phase and in the normal phase are summarized in Table III. 

The thermodynamical critical field is determined by the condition (ri) n = (ri) s and we 
obtain that [46] 



6a 2 



Hc ~ 2 V h + 3b 2 



esc 9 



(6.30) 



A first order phase transition from the super phase to the normal phase occurs for a type I 
superconductor when the external field exceeds this magnitude. 

To determine the type of a color superconductor, we repeat the analysis in [44] for a 
metallic superconductor by considering a domain wall, i.e. a planar interface between the 
super phase with CFL and the normal phase. The bulk equilibrium is maintained by an 
external magnetic field tuned at the critical magnitude (6.30) The surface tension of the 
domain wall is defined by 



f l - (f On 



the area of the interface the area of the interface' (6.31) 
A positive a favors homogeneity and the superconductor is of type I. A negative a favors 
inhomogeneity and the superconductor is of type II. The critical value of the Ginzburg 
-Landau parameter is determined by the condition a = [47]. 

On writing $ = $ + <&iT l , we observe that there are no terms in r\ that are linear 
in <3>i,...$7 or Ai,...A 7 . A consistent solution to the Ginzburg-Landau equation that im- 
plements the maximum symmetry maintained by the external conditions amounts to set 
$i = ... = $7 = and Ai = ...A 7 = 0. Taking the coordinate system with the external 
magnetic field in z-direction and the interface parallel to yz-plane, the solution ansatz is 
taken to be 



1 2 

$ = g^CFL^ + 2v) + "y|0CFL(M - v)f 8 



(6.32) 



V 



^—3a 
9 



Ay cos 6 
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with u, v and A functions of x only. The boundary condition reads 



{J [° r ^°° , (6-33) 

10 tor x — > — oo 

and 

fO for x — ► oo Cfi 

I —1 for x — > — oo 

In terms of the dimensionless field variables u, v and A, and the dimensionless coordinate 
s — |, the surface tension becomes 



6a_ 2 
6 



A f° rfs fl(A' - l) 2 + -^(u' 2 + 2v' 2 ) + -A\2u 2 + v 2 ) -\{u 2 + 2v 2 ) (6.35) 



+ 1(2 M 4 + 2u 2 v 2 + 5v 4 ) + ±-p(u 2 - v 2 ) 2 

lo lo 



with f = where « is the Ginzburg-Landau parameter defined in (6.26) and 

b\ — 3bo 

is another dimensionless parameter and is equal to — \ for perturbative coefficients (6.10). 
The minimization of (6.35) generates three coupled equations of motion, whose solution, 
when substituted back to (6.35) gives rise to the minimum surface tension, o" m m.( K >P)- The 
critical GL parameter that separate type I from type II is defined by the condition that 

0- m in.(«,p) = 0. (6.37) 

The partial derivatives of <7 m i n with respect to k or p come only from the explicit dependence 
of (6.35) on them. The contribution from the implicit dependence on k or p through u, v 
and A drops out because of the equations of motion. Therefore we have 

f^IL) < (6.38) 

v OK > P 

and 

<da mi 



dp 



> 0. (6.39) 



Consider a special field configuration u = v , which satisfies the boundary conditions. This 
trial field configuration makes the surface tension (6.35) coincide with that of an metallic 
superconductor, for which the critical GL parameter is k = -j= ~ 0.707 [44]. But the 
equations of motion may not be satisfied with u = v everywhere and therefore o" min . < at 
this k. It follows then from (6.38) that the critical GL parameter for a color superconductor 

K c (p) < -±=. (6.40) 
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Furthermore we have 

> 0. (6.41) 



d^c = /dK\ { dp ) K 

dp V(5p/cr min .=0 ^ ^^nin. ^ 



The color-flavor locked condensate corresponds to u — v and the difference u — v rep- 
resents one component of the octet under a simultaneous color-flavor rotation. Numerical 
solution to the equations of motion shows that this octet component does show up near the 
interface while the CFL condensate occupies the bulk x > 0. For weak coupling, (p = — |) 
we find that [46] 

k c ~ 0.589. (6.42) 

In terms of the Ginzburg-Landau coefficients (6.10) for weak coupling, we have a type I 
color superconductor if 

k B T c < (6.43) 

2 

and a type II color superconductor otherwise, where as = f^- 

D. Derivation of the Ginzburg-Landau theory from QCD one-gluon exchange. 

In this subsection, we shall sketch the main steps that lead to the Ginzburg-Landau free 
energy functional (6.4) with the coefficients given by (6.10) [29]. 

Starting from the path integral representation of the free energy of a quark matter, 

exp ( - (3F) = J [dA] [dijj] [diP}e- SE , (6.44) 

where 



IE 



o 



dr / (Tr£ (6.45) 



with C the QCD Lagrangian (2.1). The color superconducting order parameter can be 
triggered by the external source 

AS = J2 'Tr[J|(P)Of (P) + Ol(P)J,:(P)], (6.46) 

k,P,h 

where 0%(P) and O^(P) are the quark bilinear forms of helicity h, i.e. 

Pi(P)]?£=a%*(-P-)<%*(P+) (6-47) 

and 

mmi = * c L h (p + wL h (-p-)- (6.48) 
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with the color-flavor indices suppressed in (6.46) and the trace Tr is defined the same way 
as for \1/ in the beginning of this section. The capital letter P ± denotes the four-momentum 
(±|+p, — v). The Grassmann numbers a and a are related to the quark field in the coordinate 
representation, ip and ip through 



a c f , h ( p ) = V ^ jf ^e-^ + -4(p>^(r, r) (6.49) 

with Uh{p) the positive energy solution to the Dirac equation of helicity h. The summation 
over the four momentum P in (6.46) extends only half p-space. The pairing involving 
antiquarks is neglected. 

Upon substitution of S E + AS for S E in (6.44), we have F —> F + AF with 

exnf BAF) MMft^ (a 50 ) 

exp(-/3AF) - m]mme s E ■ ( 6 - 5 °) 

Both AF and the induced order parameter 

k[ >~ J[dA][d^[di>}e-s E ^ bL) 

and B* K (P) can be expanded according to the powers of J and J*. The Ginzburg-Landau free 
energy functional is obtained by the power series expansion of the Legendre transformation 
of AF, 

T = AF - - £ 'Tr[ J|(P)P fe -(P) + B* S (P)J S (P)] (6.52) 
P k,P 

up to the forth power of B%(P). Retaining only the color anti-triplet component of B^(P), 
we find that 



r = - £ < P'\M U \P > Mt(P)S £ (P) + ^ £[^ 2 + (P " ^) 2 ]Tr[4(P)%(P)] 2 , (6.53) 
where 

< P'\M^\P >= -5- 1 (P + )5- 1 ( - P-)^, P + f ,(P'|P) (6.54) 

v 2 

with S(P) the full quark propagator and rg(P'|P) the 2PI vertex function for diquark 

2 

scattering, all referring to positive energy states. The dependence of the quartic coefficient 
on the interactions and on the total momentum is neglected. 

The kernel Aig is isomorphic, upon multiplying the product of quark propagators on 
the right, to 1 — K^s with the kernel of the Dyson-Schwinger equation, eq.(5.10) at 
q — | and 5 = 0. Introduce the normalized eigenfunction u%(P) of Ai^ with the minimum 
eigenvalue £, i.e. 
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4^E < P \M U \P' > u u (P>) = £u n (P), (6.55) 



we have 

8tt 2 /m I., 



. k 2 B T c (T- V;,) • -/.-'• (6.56) 

with T c the same critical temperature determined before. Since the Ginzburg-Landau theory 
is expected to be valid when |T — T c \ << T c and the two terms on RHS of (6.56) becomes 
comparable, we have k « ksT c and the angular momentum mixing discussed in the context 
of LOFF pairing, which occurs when k ~ k B T c , can be neglected. 
On writing 

B S (P) = V6*MP) (6.57) 

with 

m^ = -^=J d 3 re- i%r y(r) (6.58) 

and substituting (6.57) to (6.53) we obtain the Ginzburg-Landau free energy functional (6.4) 
with the coefficients (6.10). The form of the gauge coupling is dictated by the requirement 
of gauge invariance. 

The same expression of the Ginzburg-Landau coefficients in terms of fi and T would be 
obtained if the underlying dynamics is given by a NJL effective action. 



VII. CONCLUDING REMARKS. 

In this lecture, I have reviewed the perturbative aspect of the color superconductivity 
in the regime of asymptotic freedom. The approximations made so far are systematic and 
the result obtained are theoretically important. But there is not a known mechanism which 
can maintain a quark matter with a baryon density such that ji » Aq CD at equilibrium 
against gravitational collapse. In the core of a neutron star, the baryon density is not likely 
to exceed that of a normal nuclear matter by one order of magnitude and the corresponding 
chemical potential is speculated to be few hundreds of MeV. Taking u = 400MeV and 

2 

Aqcd = 200MeV as a benchmark, we have as = — I- It follows from (2.8) for iV c = 
Nf = 3 and J = that T c = 3.5MeV. The Ginzburg-Landau parameter falls in the type I 
region according to (6.43). But this estimations may subject significant corrections for the 
following reasons: 
1. The higher order corrections: 

The sub-sub-leading term of the perturbation series (3.1) takes the form X'g with the 
ansatz 

X' = c\ng + c. (7.1) 
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The first term remains infrared, but the second term entails a matching between the infrared 
contributions to ultraviolet ones and the result depends on the definition of g. It is important 
to notice that the diagrams which seems contributing according to explicit powers of g 
turns out not so if there is no forward singularity. The crossed box diagram discussed in 
subsection IV.E is an example. While it will take major efforts to collect all contributions to 
the coefficients c and d, some of them can be obtained readily and their magnitudes serve 
an clue to the accuracy of the perturbative expansion when extrapolated to the realistic 
baryon density. One contribution to c comes from the coupling constant g under the non 
Fermi liquid logarithm of the quark-self energy. In terms of the scaled Matsubara energy £>, 
the logarithm of AK in the subsection IV. B becomes 



Following the same steps there, the second term contributes to c a term A contribu- 
tion to d comes from the quasi-particle damping rate and has been calculated in [18]. It 
contributes to d a term equal to Including these two corrections only, the transition 

temperature may be written as 



The d term of (7.1) may also be absorbed into a redefinition of Aqcd- 

The instanton mediated pairing force considered in subsection II. A, which is suppressed 
in the region of asymptotic freedom may come to play significant roles for the realistic 
baryon density. Recently the author of [48] found that the the gap energy produced by the 
one-gluon exchange and the instanton induced pairing force is considerably enhanced for a 
moderate chemical potential relative to that given by the one-gluon exchange alone. 
2. Exotic pairing state 

At a intermediate chemical potential, say fi = 400MeV the quark masses and their 
difference among different flavors has can no longer be ignored. Because of the large mass of 
s quark and the restriction of charge neutrality, we are in the situation of two flavor pairing 
in a mismatched Fermi sea. One candidate is the LOFF pairing discussed in section V, which 
gives rise to a crystalline structure of the long range order and may shed new lights on the 
glitching phenomena of a neutron star [38] . The other possibility is spin-one pairing between 
two quarks of the same flavor, which does not suffer from the mismatch [24] [25] [26]. The 
energy scale of spin-one pairing is, however very low, of the order of keV and it may occur 
in the late stage of the cooling history of a neutron star [21]. Some of the spin-one pairing 
states exhibit an electromagnetic Meissner effect [25]. 

A potential competing pairing state is the gapless superconductivity suggested recently 
[49] [50] [51] [52]. In the presence of a Fermi momentum difference 5 (defined in section V ) 
of pairing quarks, the solution to the gap equation with zero net momentum, plotted against 



In — + 3g + const.. 



(7.2) 




(7.3) 



for N c = 3, where 




(7.4) 
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8, contains two branches. One of them, with A > 8, gives rise to gapped excitation spectrum 
as usual and the other branch, with A < 8, gives rise to gapless spectrum. Regarding the 
free energy as a function of the gap parameter A, and 8, the solutions along the second 
branches are saddle points and the corresponding state is unstable without constraints. But 
as shown in [51] for NJL effective action with moderate pairing strength, the constraint of 
charge neutrality, however corresponds to a trajectory on A — 8 plane that intersects the 
gapless branch only. The free energy is minimized at the gapless solution of the gap equation 
along the trajectory of charge neutrality. The thermodynamics of this pairing state is largely 
characterized by the gapless feature of the excitations. 

Other candidate color superconductivity in the presence of a mismatched Fermi sea such 
as the gapless state with CFL [53] and the mixed state of normal and super phases [54] have 
also been considered in the literature. 

Color superconductivity continues to be an active field of research. There are many 
theoretical issues still to be addressed. The experimental confirmation its existence depends 
on the progresses along the following three avenues: 1) A robust extrapolation of the first 
principle calculations to moderate chemical potential by incorporating nonperturbative ef- 
fects. 2) A practical simulation method that can handle the fermion sign problem. 3) A list 
of clear cut signals to identify CSC phase in a compact steller objects. With joint efforts 
of high energy physicists, nuclear physicists and astrophysicists, it is quite feasible that a 
much deeper insight can be gained on the high density area of the QCD phase diagram in 
near future. 
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APPENDIX A 

In this appendix, we shall derive the formula for the group theoretic factor of the scat- 
tering amplitude of two quarks via one-gluon exchange process specific to an irreducible 
representation of di-quark system. Each of them can be in any irreducible representation of 
the gauge group. 
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Denoting the representation space of the two interacting quarks by R\ and R2, the prod- 
uct of which can be decomposed into a number of irreducible representations, say i?'s. The 
group theoretic factor associated to the one-gluon exchange diagram specific to a particular 
di-quark representation, R is an eigenvalue of the operator T£T£ which acts in the product 
space R\ <S> R2, where 

n = t c Ri ® 1 (ai) 

and 

It = 1 ® T C R2 (A2) 

with Tft the generator of the representation R. The negative eigenvalue gives rise to pairing 
while the positive one to repulsion. As only the quadratic form of the generators is involved, 
we may follow the same methodology of diagonalizing the scalar product of two angular 
momentum of SU (2) group, which amounts to 

Ji-J 2 = \{P ~ J? - Jt) = \W + 1) - Wi + 1) - W2 + 1)] (A3) 

— * — * — * 

with J = Ji + J 2 - 

The group generator in the product space R\ ® R2 is 

T c = rpc + T c 

A decomposition parallel to (A3) reads 

T°T% = - [T C T C - T{T{ - T 2 C T 2 C ] ( A5) 

and the eigenvalue of the operator TfTlf specific to the diquark representation R is 

\{C R -C Rl -C R2 ) (A6) 
with C R the second Casmir in the irreducible representation R. 

APPENDIX B 

In this appendix, we shall evaluate explicitly the Fredholm determinant pertaining to 
the integral equation (2.9), i.e. 

00 

v = E(-) nfc2n ^ ( B1 ) 


where 

1 r b r b 

T> n = — dx n ... dx 1 D n (x 1 , ...,x n ) (B2) 

n\ J a J a 
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with D n (xi, ...,x n ) the determinant 



D n (xi, x n ) 



X\ min(xi, X2) min(xi, X3) 

min(x2,Xi) X2 mm(x2,X3) 

min(x 3 ,Xi) min(x 3 ,a;2) x 3 



min(x n ,xi) min(:r ra ,X2) min(a; n ,j; 3 ) 
In terms of the variables 



mm{xi,x n ) 
min(x 2 ,x n ) 
min(x 3 ,a; n ) 



(B3) 



b — a 



we have 



where 



(b-a) n I n + a(b-a) n - L J n , 



/. = !/' 

n! Jo 



and 



1 f 1 f 1 

J n = l TTT / *;„.../ ^1-^(6, f n ). 

(n — Ij! jo jo 



(B4) 



(B5) 



(B6) 



(B7) 



with -D*(£i, £ n ) the determinant of the matrix obtained from D n (£i, £„) upon replacing 
its elements of the first column by ones, i.e. 



1 min(fi,f 2 ) min(fi,£ 3 ) 

1 £ 2 min(^ 2 ,^3) 

1 min(£ 3 ,£ 2 ) 6 

1 min(£ n ,£ 2 ) min(£ n ,£ 3 ) 



min(£ 2 , £ n ) 
min(£ 3 , f n ) 



(B8) 



The permutation symmetry of the integral with respect to xi,...x n is employed in deriving 
(B5) from (B2). 

To evaluate I n , we order the integration variables such that £1 < £ 2 < ••• < Cn- We find 
that 



(B9) 



Multiplying the second column by g- and subtract the result from the first column, we obtain 
the relation 





6 


6 


6 • 


• 6 






6 






Ai(£l, £ra) — 




6 








6 


6 


6 • 





A»(£L,6»-»Zn) = ^(l - ^)D„_i(6,-,Cn) 

X ?2 y 



(BIO) 
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Recursively, we end up with 

n n—1 



A,(ei,-,&) = n&n&+i-&)- ( Bn ) 

3=1 3=1 



Introduce a new set of integration variables 



6 = VlV2---Vn (B12) 
6 = V2V3---Vn 



we have 

„\ n n—1 -i 

To evaluate J n , we order the integration variables such that £ 2 < ••• < in- Following the 
same treatment of D n , we find -D n (£i, £2, £n) = if £1 > £2 and 

D n (^-U) = Dni ^^ n) (B14) 
si 

if £1 < £ 2 . In terms of the new variables (B12), we end up with 

Jn = 7 —^7- (B15) 

(2n-l)! V ; 

Substituting (B13), (B15) and (B5) to (Bl) we obtain that 

V = cos k(b — a) — ka sin k(b — a), (B16) 
which gives rise to the same eigenvalue condition as (3.43). 

APPENDIX C 

In this appendix, we shall discuss the color sextet component of the condensate. Upon 
decomposing the diquark condensate into color 3 component 0, and 6 component Xi 

^ = <P + X (CI) 

the quadratic term of the Ginzburg-Landau free energy takes the form 

r 2 = Tr(a0t0 + aV X ) (C2) 

with a = a(T — T c ) and a' > 0. The quartic term T 4 contains all possible combinations of 
and x that is invariant under the symmetry group (2.19). Among these combinations, 
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there exits a term of the form x*0*00 and its complex conjugate. Indeed, it follows from 
the decomposition rules of SU (3) group 

3 <g> 3 = 3 © 6 (C3) 
3 © 3 = 3 © 6 
6 © 3 = 3 © 15 

that the product representation of 0*00 consists of one 6 which forms an invariant with x- 
Since this term is linear in x, a nonzero expectation value of x will be induced by a nonzero 
value of upon minimizing the Ginzburg-Landau free energy at T < T c and its contribution 
to the free energy in the super phase is smaller than that of the anti-triplet by one power 
of Tc ~ T . For the Ginzburg-Landau coefficient in the weak coupling, (6.10), we find that the 
sextet component of the condensate corrects the bulk free energy density by the amount 

Ar mir , _ 558C(5) u 2 Pl *(Tc-Ty 



50 



REFERENCES 



[1] B. Barrois, Nucl. Phys. B129 390 (1977). 

[2] S. Frautschi, Proceedings of the Workshop on Hadronic Matter at Extreme Energy Den- 
sity, N. Cabibbo, ed., Erice, Italy (1978). 

[3] D. Bailin and A. Love, Phys. Rep. 107 325 (1984). 

[4] M. Alford, K. Rajagopal and F. Wilczek, Phys. Lett. B422 247 (1998). 

[5] R. Rapp, T. Schaefer, E. V. Shuryak, M. Velkovsky, Phys. Rev. Lett. 81 53 (1998). 

[6] M. Alford, K. Rajagopal and F. Wilczek, Nucl. Phys. B537, 443 (1999). 

[7] M. Gyulassy, nucl-th/0403032. 

[8] D. T. Son, Phys. Rev. 59, 094019 (1999). 

[9] T. Schaefer and F. Wilczek, Phys. Rev. D60, 114033 (1999). 
[10] R. Pisarski and D. Rischke, Phys. Rev. Lett. 83 37 (1999). 
[11] R. Pisarski and D. Rischke, Phys. Rev. D61 051501 (2000). 
[12] D. K. Hong, Phys. Lett. B473 118 (2000). 

[13] D. K. Hong, V. A. Miransky, I. A. Shovkovy and C. R. Wijewardhana, Phys. Rev. D61, 
056001 (2000). 

[14] W. Brown, J. T. Liu and H. C. Ren, Phys. Rev. D61, 114012 (2000). 

[15] S. D. H. Hsu and M. Schwetz, Nucl. Phys. B572, 211 (2000). 

[16] W. Brown, J. T. Liu and H. C. Ren, Phys. Rev. D62, 054016 (2000). 

[17] W. Brown, J. T. Liu and H. C. Ren, Phys. Rev. D62, 054013 (2000). 

[18] C. Manuel, Phys. Rev. D62 114008, (2000). 

[19] Q. Wang and D. Rischke, Phys. Rev. D65 054005 (2002). 

[20] De-fu Hou, Qun Wang and D. Rischke, hep-ph/0401152. 

[21] K. Rajagopal and F. Wilczek, hep-ph/0011333, in At the Frontier of Particle 
Physics /Handbook of QCD, B. L. Ioffe Festschrift, edited by M. Shifman (World Sci- 
entic, Singapore); M. Alford, Ann. Rev. Nucl. Part. Sci., 51, 131 (2001); T. Schaefer, 
hep-ph/0304281; D. Rischke, Prog Part. Nucl. Phys. 52 197(2004); R. Casalbuoni and 
G. Nardulli, hep-ph/0305069. 

[22] 't Hooft, Phys. Rev. D14, 3432 (1976). 

[23] M. Shifman, A. J. Vainshtein and V. I. Zakharnov, Nucl. Phys. B163, 46 (1980); E. 

Shuryak, Nucl. Phys. B203, 140 (1982); D. I. Diakonov and V. Y. Petrov, Nucl. Phys. 

B245, 259 (1984). 
[24] T. Schaefer, Phys. Rev. D62, 094007 (2000). 

[25] A. Schmitt, Qun Wang and D. Rischke, Phys. Rev. Lett., 91, 242301 (2003). 
[26] M. G. Alford, J. A. Bowers, J. M. Cheyne and G. A. Cowan, Phys. Rev. D67, 054018 
(2003). 

[27] N. Evans, J. Hormuzdiar, S. D. H. Hsu and M. Schwetz, Nucl. Phys. B581, 391 (2000). 
[28] T. Schaefer, Nucl. Phys. B575, 269 (2000); I. Shovkovy and L. C. Wijewardhana, Phys. 

Lett. 470B, 189 (1999). 
[29] I. Giannakis and H. C. Ren, Phys. Rev. D65 054017 (2002). 
[30] M. Alford, J. Berges and K. Rajagopal, Nucl. Phys. B571 (2000). 
[31] K. Rajagopal and F. Wilczek, Phys. Rev. Lett. 86, 3492 (2001). 

[32] D. T. Son and M. A. Stephannov, Phys. Rev. D61, 074012 (2000); D62, 059902 (2000). 
[33] C. Manuel and M. H. G. Tytgat, Phys. Lett. B479, 190 (2000) 



51 



[34] M. Le Bellac, Thermal Field Theory (Cambridge University Press, Cambridge, England, 
1996). 

[35] E. T. Whittaker and G. N. Watson, A Course of Modern Analysis, 4th ed. (Cambridge 

University Press, Cambridge, England, 1962). 
[36] T. Holstein, R. E. Norton and P. Pincus, Phys. Rev. B6 2649(1973); M. Yu. Reizer, 

Phys. Rev. B40, 11571(1988), Phys. Rev. B44, 5476(1991). 
[37] T. Schaefer, Nucl. Phys. A728 251(2003). 

[38] J. A. Bowers, J. Kundu, K. Rajagopal and E. Shuster, Phys. Rev. D64, 014024 (2001). 
[39] A. J. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Teor. Fiz 47 1136(1964) (Soviet Physics 

JETP 20 762(1965); P. Fulde and R. A. Ferrell, Phys. Rev. 135 A550(1964). 
[40] S. Takada and T. Izuyama, Prog. Theor. Phys, 41 635(1969). 
[41] A. K. Leibovich, K. Rajagopal and E. Shuster, Phys. Rev. D64, 094005 (2001). 
[42] I. Giannakis, J. T. Liu and H. C. Ren, Phys. Rev. D66 031501 (2002). 
[43] J. Bowers and K. Rajagopal, Phys. Rev. D66 065002 (2002). 
[44] V. L. Ginzburg and L. D. Landau, Zh. Eksp. Teor. Fiz. 20 1064 (1950). 
[45] K. Iida and G. Baym, Phys. Rev. D63 074018 (2001); Phys. Rev. D66 014015 (2002). 
[46] I. Giannakis and H. C. Ren, Nucl. Phys. B669 462 (2003). 
[47] A. A. Abrikosov, Soviet Physics JETP 5 1174 (1957). 
[48] T. Schaefer, hep-ph/0402032. 

[49] W. V. Liu and F. Wilczek, Phys. Rev. Lett. 90 047002 (2003). 
[50] Mei Huang, Pengfei Zhuang and Weiqin Chao, Phys. Rev. D67, 065051 (2003). 
[51] I. Shovkovy and Mei Huang, Phys. Lett. B564, 205 (2003); Nucl. Phys. A729, 835 
(2003). 

[52] A. Mishra and H. Mishra, Phys. Rev. D69, 014014 (2004). 
[53] M. Alford, C. Kouvaris and K. Rajagopal, hep-ph/0311286. 

[54] P. F. Bedaque, H. Caldas and G. Rupak, Phys. Rev. Lett. 91, 247002 (2003); H. Caldas, 
hep-ph/0312275. 



52 



